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INTRODUCTION 

In  19U1,  Borrmann  observed  anomalously  large  diffracted  and  transmitted 
x-ray  beams  obtained  from  propagation  through  thick  single  crystals. ^^^* 
This  effect  has  since  been  named  the  Borrmann  effect  and  is  the  first 
published  reference  to  anomalous  transmission. 

Ordinarily  one  would  expect  x-rays  to  be  absorbed  according  to  the 
usual  expression  for  absorption:  I  =  1^jhp{-^^'*)    •     If  one  were  to  use  a 
crystal  of  such  a  thickness  that  ^Jcji-^io    ,  then  the  normally  expected 
incident  beam  would  be  attenuated  by  j^TCp  (-(o")  in  passing  through  the  crystal 
or,  in  other  words,  there  would  be  no  observable  transmitted  beam.  The 
fact  that  there  are  large  transmitted  beams  would  seem  to  indicate  that  the 
normal  photoelectric  or  kinematic  effect  is  not  valid  in  this  case.  The 
explanation  for  this  phenomena  lies  in  the  dynamical  theory  of  x-ray  diffrac- 
tion. 

Consider  figure  1.  A  single  crystal,  cut  in  the  form  of  a  parallel 
sided  slab  with  the  planes  to  be  used  in  the  diffraction  perpendicular  to 
the  slab  faces,  is  rotated  about  an  axis  perpendicular  to  the  plane  of  the 
figure.  If  the  crystal  is  oriented  in  such  a  way  that  the  Bragg  condition, 
X  =  2<^sin  ©^  is  not  satisfied,  then  the  intensity  is  given  by  I  =  I  exp(-M,f). 
If  two  crystals  are  considered,  one  for  which  p/<<i   and  the  other  having 
l^^y^'io,   then  a  rotation  of  one  of  the  crystals  through  Q^  results  in 
the  observed  phenomena  in  Part  C  of  the  figure.  For  the  thin  crystal  there 
is  a  decrease  in  transmitted  intensity  as  the  crystal  reaches  the  Bragg 
angle,  however,  for  the  thick  crystal  there  is  a  marked  increase  in 


*Ref  erences  are  on  page  5"  7 
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transmitted  intensity.  That  this  effect  is  due  to  diffraction  is  borne  out 
by  figure  2. 

If  the  thin  crystal  in  this  figure  is  set  for  Bragg  diffraction  and  a 
film  is  placed  behind  the  crystal,  there  will  appear  three  small  dark  areas 
on  the  film.  Area  (l)  is  due  to  diffraction;  area  (2)  is  approximately 
equal  in  intensity  to  area  (l)  and  at  a  distance  from  area  (l)  corresponding 
to  the  proper  angle  a©^  ^o^  diffraction;  the  third  area,  directly  in  line 
with  the  incident  beam,  is  very  weak.  The  fact  that  area  (2)  is  almost  as 
intense  as  area  (l)  and  that  the  distance  from  area  (2)  to  area  (3)  is 
proportional  to  the  thickness  of  the  crystal  would  seem  to  indicate  that  the 
x-rays  are  channeled  down  the  atomic  planes  and  are  attenuated  little  by 
the  crystal. 

The  dynamical  theory  of  diffraction  explains  this  phenomenon  well,  due 
to  the  fact  that  in  it  diffraction  is  handled  as  being  a  scattering  effect 
of  the  crystal  as  a  whole,  whereas  the  kinematic  theory  fails  to  predict 
the  result  because  it  considers  the  scattering  as  being  a  process  associated 
with  individual  atoms  in  the  crystal.  The  dynamical  treatment  may  be 
considered  as  looking  at  the  crystal  as  though  the  atom  planes  along  the 
direction  of  the  x-rays  were  mirrors.  At  the  correct  wave  length  a  wave 
field  is  established  and  little  or  no  attenuation  takes  place.  If  the 
crystal  is  too  thin,  there  are  not  enough  atomic  planes  to  establish  this 
x-ray  wave  field  so  a  transmitted  beam  is  reduced  at  the  Bragg  angle. 

This  anomalous  transmission  does  not  take  place  along  all  directions 
in  the  crystal  to  the  same  degree.  Various  parameters  have  to  be  taken 
into  consideration  to  determine  the  amount  of  anomalous  transmission  along 
a  given  direction.  This  amount  of  transmission  is  related  to  a  quantity 
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called  epsilon,  e,  which  varies  from  0  to  1  depending  on  the  direction  of 
the  x-rays  in  the  crystal.  The  quantity  e  may  be  looked  at  as  being  a  term 
added  to  the  ordinary  absorption  such  that  1  =  1  exp  l_-u.t{i-e)]  .     That  is, 
if  e  =1,  then  I  =  I^,  and  if  e  =  0  the  standard  absorption  results. 

This  quantity  epsilon  is  not  only  related  to  the  direction  in  the  crystal, 
but  may  also  be  related  to  temperature  effects,  i.e.,  lattice  vibration. 
B.  Okkerse-'--'-  has  shown  that  the  temperature  dependence  of  the  quantity 
epsilon  is  e  =  ZQey.p(rf^)     .  Epsilon  is  a  measured  value;  Eq  is  the  calculated 
value  related  to  the  distribution  of  absorbing  power  in  the  crystal  of 
the  various  atoms  as  obtained  from  the  incoherent  scattering  factors  of 
the  electron;  M  is  the  Debye-Waller  factor  M  =  ^   <*-<.'■>  -iin^    (  u. 
is  the  root  mean  square  amplitude  of  vibration). 

It  is  the  purpose  of  this  experiment  to  determine  the  values  of  epsilons 
for  several  calcite  (CaCO^)  planes  both  theoretically  and  experimentally  and 
make  a  comparison  between  them. 


THEORETICAL  CALCULATIONS 

The  theoretical  values  for  epsilon,  e^,  can  be  calculated  in  two 
different  ways.  Both  methods  will  be  discussed  here. 

The  first  method  of  calculating  epsilon  is  based  on  a  calculation 
made  by  Boris  W.  Batterman-^  on  germanium. 

Nomenclature 

F(hkl)  =  F' (hkl)  +  ZF"  (hkl)  Crystal  structure  factor 

P    Polarization  factor  =  1  or  cos  29  for 
<^  and  n  polarization,  respectively 

Og       Bragg  angle 

H,  j^   Reciprocal  lattice  vector  hkl  reflection 

f  .,f .  Atomic  scattering  factor  for  the  j"^^  electron 

■^  in  the  atom  with  and  without  dispersion,  respectively 

2?  Number  of  electrons  of  type  j  in  atom 

7.1, '^i   Wave  length  and  angular  frequency  of  absorption 
edge  for  j  electrons 

Tq       Classical  electron  radius 

N   Number  of  unit  cells  per  unit  volume 

X   X-ray  wave  length 

M   Debye  Waller  factor  for  amplitude 

As  was  stated  in  the  introduction,  the  dynamical-  theory,  unlike  the 

kinematic  theory,  takes  into  consideration  scattering  of  x-rays  by  the  entire 

crystal  rather  than  by  point  sources.  Consequently,  in  the  dynamical  theory 

Maxwell's  equations  are  solved  for  the  condition  of  a  complex  dielectric 

constant  within  the  crystal.  A  resultant  equation  for  the  ratio  of  the 

diffracted  to  the  incident  intensity  is  derived  with  some  calculation.  No 

attempt  will  be  made  to  derive  this  equation.  A  good  explanation  of  this 


derivation  may  be  found  in  Zachariasen,-^'^  or  Batterman  and  Cole. 

For  the  case  of  centrosymmetric  crystal,  this  value  for  the  ratio 
of  the  intensities  is  given  by  Hirsch"  as 

-£a  =  -kZ-    :       1  jy,,  (-  He  t,/^  \   t^^jL  a.  A  K 

where  H^    is  the  cos  of  the  angle  of  incidence  to  the  normal  of  the  incident 
surface,  /y,  is  the  linear  absorption  coefficient,  and 

From  this  equation  a  "new"  absorption  coefficient  can  be  derived. 

It  is  necessary  that  an  average  over  the  two  types  of  polarization  be 
taken  letting 

or  equating  it  to 

results  in 

epsilon,  e,  is  defined  to  be  the  ratio 

e  ^  Pf"(A'»^3/p"(°°o) 
where  this  is  the  ratio  of  the  imaginary  parts  of  the  scattering  factor  for 

the  hkl  and  000  reflections,  respectively.  Thus,  the  theoretical  value  of 

epsilon  can  be  obtained  by  finding  F"  (hkl)  and  F"  (OOO),  respectively. 

The  structure  factor  for  a  crystalline  (hkl)  reflection  is  given  by 

the  summation  over  all  atoms  in  the  unit  cell 

where  r  is  a  vector  from  an  origin  in  the  unit  cell  to  the  n"*^^  atom.  All 
possible  values  of  a^p  (air/  H.    •-^.«)  are  calculated  for  each  atom  in  the 
crystal.  The  position  of  the  atoms  is  obtained  from  Wycoff  .•'•3 
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The  atomic  scattering  factor  /  for  each  atom  may  be  considered  to  be 

imaginary,  that  is 

f  =  fO  +  Af  +XAf" 

where  Af  and  Af"  are  the  Honl'^  or  dispersion  corrections  and  are  dependent 
on  the  frequency  of  the  incident  radiation,  while  f°  is  the  form  factor 
depending  on  charge  distribution  and  is  independent  of  frequency.  For  a 
particular  crystal  f°,  Af,  and  Af"  vary  with  a  change  in  ^^  g:  ;  that  is, 
if  two  planes  of  diffraction  have  different  d-spacing,  the  values  of  f°, 
Af,  and  Af"  are  not  the  same  for  the  two  planes.  The  term  Af  is  the 
modification  of  the  atomic  scattering  when  the  radiation  used  is  close  to 
the  absorption  edge  of  one  of  the  atoms  in  the  sample.  Since  the  wavelengths 
that  were  used  here  are  not  near  the  absorption  edges  for  any  of  the  atoms 
in  calcite,  Af  is  small  compared  to  f°  and  will  need  no  longer  be  considered. 
So  for  the  purposes  of  this  work  f  reduces  to 

f  =  f  0  +  <f  Af " 

The  value  of  f °  is  given  by 

0 

fo  =  ^^   ?^.  i^ 

Where  ?.  is  the  number  of  electrons  of  a  given  type  in  the  atom  and  f ^  is 
the  scattering  factor  for  the  electron.  This  scattering  factor  is  derived 
from  the  Thomson  classical  scattering  of  x-rays  by  the  electron.-^ 

The  x-rays  incident  upon  the  atom  give  rise  to  an  oscillating  electronic 
dipole  movement  which  radiates  x-rays.  To  get  the  total  amplitude  scattered 
by  the  j"*^^  electron,  the  dipole  radiation  scattered  from  the  fraction  of  the 
2^^   electron  in  the  volume  increment  dv  is  s\immed  over  all  space  taking  into 
account  the  path  length  differences  between  the  elements  dv.  Therefore, 


u 


where  <f>   is  the  phase  angle  depending  upon  the  angle  of  scattering,  f   is  the 
number  of  electrons  in  the  element  dv.  D  is  the  amplitude  scattered  by 
that  charge.  D-  is  dependent  on  electron  type  only  so  it  may  be  taken  out 
of  the  integral.  This  gives 

U-   ^'Jf  where   /J-  j"  ^  ^      ■ 

Parrattsnd  Hempsteadl2  showed  that 

o.  -  0  -^  /,  +  -  ^4 )  ^  (  ^-f.  ^"'^  ^  ,  . 

where  co  is  the  angular  frequency  of  a  virtual  oscillation  associated  with 

the  y^^   electron,  for  which  there  are  W.(a))da)  oscillators  in  the  frequency 

range  dco,  Mq  is  the  incident  x-ray  frequency,  w  is  the  frequency  of  the 

j'tn  electron  absorption  edge,  and  k^  is  a  damping  constant.  Combining  the 

two  expressions  for  f •  and  the  expression  for  D-  then 

J  J  • 

Next  /Ij  must  be  found  for  each  electron.  Again  from  Parratt  and  Hempstead' s^^ 
treatment  of  anomalous  dispersion  it  is  found  that  for  Cu  K"^    radiation 

where  g.  is  the  oscillator  strength  for  all  the  3  electrons,  X.  is  the 
absorption  edge  wavelength  and  P^  =  ll/U,  2,  5/2,  5/2  for  the  K,  L-^,  Lg, 
L^  and  M  absorption  edges,  respectively.  The  values  of  f .  for  the  various 
electrons  in  calcium,  carbon,  and  oxygen  v^ere  obtained  from  three  papers 
by  Freeman, 15  (Ca,  0)  and  Keating  and  Vineyard^S  (C) .  The  values  for  f^, 
^L-|_*  ^L2  3»  ^"^  %  ^^®  plotted  for  several  values  of  ^in  OJ^.  ,  see  figures  3, 
U,  and  5«  The  theoretical  calculations  of  epsilon  for  the  211  plane  is  used 
to  illustrate  the  procedure. 

As  will  be  shown  later,  the  crystal  structure  factor  for  the  211  plane 
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of  calcite  is  given  by 

'  0.11    "^K^c^      'c   '•  /         ^      x^  6.0^ 

The  atomic  scattering  factors  are  written  in  terms  of  their  real  and  complex 
parts  so  F211  =   a  [Ui  ^  &  a^'  (-  C  -  ^l'  -^  ^  I')  ] 


The  values  for  these  parameters  are  found  from  table  1  and  figures  3»  U>  snci 

Af^"  =  1.091  X  0.998  +  0.1h7l  X  0.968  +  0.859  X  .97U  +  0.0006  X  0.7 
Afc"  =  O.OIOU  X  0.982  +  0.0012  X  0.$79  +  0.0000^3  X  0.UU6 
Afg"  =  0.0102  X  0.990  +  0.0108  X  0.7h6  +  0.0003U  X  0.630 
^''f^a  =  1.317        A"f^  =  0.01093  A"f^  =  0.01835 


thus,  F21X  =  2.692 


-  ^211" 

^211  - iT 
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000 


From  table  2 
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TABLE  2:   CHART  OF  THE  CRYSTAL  STRUCTURE  FACTORS 


Fy^y^-^  h  +  k  +  1     mors     h     p     1     q     ^  -^  1 

2(fca  ^  ^c  ^  ^^o) 

2  (fca  ^  f  c  -  f  o) 

2(fca  -^  ^c  -  ^o) 

2(fca  ^  fc  -  ^o) 

2  Cf  ca  ^  f  c  ^  ^o) 

2(fca^fc-^o) 
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TABLE  2:       CHART  OF  THE  CRYSTAL  STRUCTURE  FACTORS   (continued) 
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TABLE  2:   CHART  OF  THE  CRYSTAL  STRUCTURE  FACTORS  (continued) 


Fyj^^l  h  +  k  +  1       m      h     p    1     q    ^^-^ 


2  (f  ca  -  f  c  "  ^c)         ^  0      0         ° 

2  (f  ca  +  f  c  -  ^o)         ^  ^      ^     °     E     0 
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h'  „^  =  2(f   •••  f   +  -^f  ) 


000 

F    =2 

000 


n,  ?^/OoU'^'.  ^*'^^'^- 


Since  for  (hkl)  =  (OOO),  .5ilL£  =  0,  then  all  the  incoherent  scattering 
factors  are  1  so 

^ooo"  "  2  [1.0910  +  O.lhTl  +  0.0359  +  0.0006  +  O.OlOh 

+  0.0012  +  0.0001  -f  3(0.0102  +  0.0108  +  0.0003)]  =  2.800 

Thus, 

F"2ii   2.692 

eo-iT  =  =  — =  0.961 

'""   F",,,   2.800 

Table  3  has  the  other  values  of  epsilon  calculated  for  this  crystal 
using  this  method. 

TABLE  3 
CALCULATED  VALUES  OF  EPSILON 


60211  0.961 

^0310  0-919 

^ou22  0-919 

eolOl  0.931 

^0213  0.93U 

^oU23  0.019^ 


-  ^'j--' ^ff-ir^-^-  ■•-- 


2h 


The  previous  calculation  of  epsilon  was  based  on  the  solution  of 
Maxwell's  equations  for  the  x-ray  case.  Epsilon  can  be  calculated  by  a 
different  procedure  which  is  based  on  the  absorption  coefficient  of  each 
individual  atom  in  the  crystal.  If  the  total  electron  field  in  the  crystal 
is  considered  as  being  the  vector  sum  of  the  incident  and  diffracted  beams, 
then 

lef2  =  E^  2  +  Eh  2  +  2E^  •  E^  exp  (-2TiiH  •  r) 

Referring  to  Figure  6,  let  r^  be  the  position  of  the  nucleus  of  an  atom 
in  the  wave  field  and  r  the  vector  to  the  volume  dv.  Let  p'(r^)dv  be  the 
contribution  of  the  element  dv  to  the  average  absorption  coefficient  of  the 
atom.  Let   /^.  =  j  ^  (r^)dv  be  the  average  absorption  coefficient  of  the 
atom.  If  it  is  assumed  that  the  absorption  coefficient  at  any  point  is 
proportional  to  |£|  ,  the  absorption  coefficient  of  the  atom  at  r  in  the 
field  can  be  written  as 

=  t[(Eo2  +  E^^)f^^{r^)d■v  =   2E^  •  Ej^J'exp(-2ni  H.r)p(r^)dv 
=  t[(Eq2  +  £^2)^^.  +  2Eq  •  Eh  exp(2ni  H-r^)  } 
XT(exp(2ni  Ht^) /^  (rg^)dv 
where  the  integral  is  over  the  volume  of  the  atom  and  T  is  a  proportionality 
constant.  If  this  value  of  /^(r^)  is  found  for  each  position  in  the  crystal, 
an  averaged  fJ.   results.  If  the  last  integral  is  abbreviated  A4m.  ,  then 


'      Eo^  ^  %2    ^°  / 


The  formal  theory  yields 
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The  similarity  makes  it  desirable  to  equate  e  and  f^  ,  and  as  it  is  defined, 
e  is  related  to  the  weighted  distribution  of  absorbing  pov;er  within  the 
atom. 

To  calculate  e  consider 


( ^\       e    -      -^  -    -L  C  iwp  (iTr/H.;«0  H^^"-')  *^ 


where  j's  are  summed  over  all  atom  positions.  As  an  example,  e  for  the 

211  plane  of  calcite  is  calculated  here. 

From  Wye of f^  the  following  positions  for  the  atoms  in  calcite  are 

obtained 

Ca  (0  0  0,  1/2  1/2  1/2)    C  (l/U  l/U  l/U,  3/U  3/U  3/U) 

0  (0  1/2  VU:  1/2, 3/U, 0;  3/U, 0,1/2; 
0,-1/2, -3/U;  -l/2,-3/U,0; 
-3/U, 0,-1/2) 

For  211,  H  •  r  for  each  atom  is  calculated  where  H  =  2b-]_  +  2b2  +  2bo  and  r's 

are  given  above 


Ca 

C 

0 

0 

1 

5/U 

-5/U 

2 

3 

7/U 

-7/U 

2 

-2 

H 


The  equation  for  epsilon  is  given  in  terms  of  absorption  coefficients 
per  atom,  whereas  available  data  are  in  terms  of  absorption  coefficients 
per  unit  mass,  so  a  conversion  from  one  set  of  units  to  the  other  is 
necessary. 


2B 


The  mass  absorption  coefficient  for  a  particular  atom  type  is 


=   f^c 


/nrr-L.^<jk. 


M, 


Since  the  absorption  coefficient  is  in  terms  of  mass  absorption  coefficients, 
these  can  be  substituted  in  the  equation  for  the  atomic  absorption 
coefficients. 


eni-s.  %  6 


A-r.   VV/:r 


1 


M, 


A^-. 


'  I    '■^^^•/^..-.(^-•'^-/) 


<«' 


As  an  example,  the  calculated  value  of  e  for  the  211  plane  using  Cullity's^ 
value  for  mass  absorption  coefficient 


=  (UP  X  172)(2)  ^   (12  X  g.$)(2)  +  (16  X  12.7)(6)  _ 


7556 


=  0.9U2 


THE  EXPERIMENT 

In  order  for  a  wave  field  to  be  established  in  the  crystal,  the  structure 
of  the  crystal  must  be  very  nearly  perfect.  If  there  is  an  imperfection, 
either  a  foreign  element  or  a  strain  producing  a  change  in  d-spacing 
different  from  the  rest  of  the  crystal,  then  the  constructive  interference  . 
necessary  for  anomalous  transmission  will  not  occur. 

Several  crystals  upon  which  similar  experiments  have  been  done  have 
been  synthetically  produced.  These  crystals  exhibit  a  high  degree  of 
perfection.  To  date  no  synthetic  crystals  of  calcite  have  been  produced 
large  enough  for  the  present  research.  Consequently,  a  naturally  occurring 
calcite  crystal  had  to  be  used  for  this  experiment.  Rather  good  calcite 
samples  were  procured  from  J.  A.  Bearden   (Johns  Hopkins  University). 

The  determination  of  the  perfection  of  a  crystal  is  made  using  the 
Lang  Topographs.  In  the  Lang  technique,  the  crystal  is  placed  on  a  camera 
that  allows  the  crystal  to  be  translated.  The  film  holder  is  also  attached 
to  the  movable  slide.  Between  the  crystal  and  the  film  there  is  a  vertical 
slit.  The  slit  is  stationary  with  respect  to  the  camera  and  x-ray  beam. 
Referring  to  the  figure  7 ,  if  the  crystal  translates  as  is  shown  and  the 
camera  is  set  at  the  proper  angle  such  that  the  Bragg  condition  is  satisfied, 
then  a  one  to  one  mapping  of  the  crystal  will  be  made  onto  the  film. 
Wherever  there  is  an  imperfection  in  the  crystal,  a  dark  area  will  appear  on 
the  film.  The  reason  that  a  dark  area  appears  wherever  there  is  an 
imperfection  is  due  to  the  break  down  of  first  order  or  primary  extinction. 
A  good  explanation  of  the  Lang  topograph  is  given  by  Lang.° 

Topographs  of  three  calcite  crystals  are  shown  in  figures  8,  9,  10,  and 
11.  The  first  few  crystals  were  quite  imperfect,  notably  the  B  and  L 
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FIGURE  7 
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crystal,  figure  8.  However,  after  a  few  trials  it  was  noticed  that  the 
imperfection  that  appeared  in  the  Argentine  crystals  followed  the  basic 
contour  of  the  surface.  Consequently,  it  was  decided  that  an  attempt 
should  be  made  to  polish  the  crystal. 

To  do  this  polishing,  a  Buehler  polishing  apparatus  was  used,  see 
figures  12,  13,  lU  and  15 •  In  this  unit  the  crystal  is  mounted  on  the  five 
holed  wheel.  Wooden  blocks  were  inserted  in  the  hole  in  order  that  the 
crystal  could  be  glued  to  it.  The  five  holed  piece  was  placed  on  the 
polishing  wheel  and  the  driving  mechanism  was  placed  over  this.  The 
driving  gear  counterrotated  the  sample  wheel  as  the  polish  progressed.  The 
Buehler  instruction  manual  calls  for  a  pressure  of  fifty  to  sixty  pounds 
to  be  placed  on  the  wheel.  However,  it  was  realized  rapidly  that  this  was 
impractical  since  cloths  were  ripped  and  crystals  cracked.  Consequently, 
all  crystals  were  polished  without  adding  to  the  pressure  developed  by  the 
weight  of  the  wheel. 

The  crystals  were  brought  down  to  the  desired  thickness  using  No.  600 
grit  paper.  The  rough  polishing  was  accomplished  first  with  Billiard  cloth 
and  6  micron  aluminum  oxide  polish,  follovred  by  micron  cloth  and  1  micron 
aluminum  oxide.  The  final  polish  was  with  0.05  micron  Gamma  Al-umina  on  a 
velvet  cloth.  The  entire  procedure  took  approximately  tvjo  hours  per  crystal 
and  produced  a  highly  polished  surface.  Some  of  the  crystals  were  etched, 
but  this  did  not  show  up  on  the  topograph,  so  the  later  crystals  were  not 
etched. 

After  the  polishing,  several  crystals  exhibited  almost  perfect 
topographs,  figure  11.  The  perfect  topograph  has  only  the  edge  of  the  crystal 
apparent,  and  the  center  of  the  crystal  is  either  clear  or,  at  least. 
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contains  a  uniform  light  density.  Only  crystals  that  exhibited  this 
condition  and  were  of  sufficient  size  to  be  handled  easily  (approximately 
1  cm^)  were  used. 

As  was  stated  earlier,  in  order  that  a  crystal  exhibit  anomalous 
transmission,  it  must  have  a  jJ,x  >>  I      that  is,  the  linear  absorption 
coefficient  times  the  thickness  must  be  much  greater  than  one.  However,  it 
should  not  be  too  large  because 

so  if  jjot    is  too  large,  say  on  the  order  of  15,  then  only  x-rays  from 
extremely  strong  diffracting  planes  with  values  of  e^  1  will  be  seen.  As 
it  was  desired  to  work  with  several  planes  of  low  t,  a   y.Jt    of  about  five  was 
desired.  With  this  in  mind,  a  0.^13  mm  crystal  was  chosen  for  calcite  in 
CuK*^  radiation  and  a  2.^8  mm  crystal  was  chosen  for  the  M^KK  radiation. 
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DETERMINATION  OF  THE  PLMES  IN  GALCITE  TO  IWESTIGATE 

Of  particular  importance  as  regards  the  planes  chosen  for  this 
experiment,  are  (1)  the  crystal  structure  factor  for  the  plane,  and  (2)  the 
angle  between  a  particular  plane  and  the  other  planes  in  the  crystal.  The 
crystal  structure  factor  is  defined  by  the  equation 

where  hkl  are  the  Miller  indices  for  the  plane,  f^^U-are  the  atomic  scattering 
factors  for  the  various  atoms  in  the  crystal,  and  u^  ,  v^  and  w^  are  the 
position  of  the  various  atoms  in  the  unit  cell  of  the  crystal. 

For  calcite  the  atoms  in  the  lattice  have  the  following  positions 
CA  C  0 

0    0    0        lA  l/U  l/U         0   1/2  3A      0   -1/2  -3/U 
1/2  1/2  1/2       VU  VU  3/U        1/2  3/U   0     -1/2  -3/U    0 

3/U   0   1/2     -3/U    0   -1/2 
Substituting  these  values  in  equation  8  for  u  ,  v     and  w 

^^KX  '  /c^L  -e/x/' («7r'i  [oA  -to^i   to^j)    4-  jL^/'fzrJQ-^/^   V-  V:^  +  "^Z*^)  J 

-»  ^^    f  Jh'P  (arr/  QjA.  ■*■  I   i^  +  Jy^  J^I  )  ■*■    •*^'"  (^'^^  ii^  *  \'^   *■  OiJ  )  J 
^^^     {J^fl(J!lTrJ[_^K  +  o/^^.uJi'])     -e  ax/»  (aff-^*  £oX  +  ^  *^  -*  i  jf]  ^  J 
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When  these  are  reduced,  the  expression  becomes 

•using  the  relation  cos  (A  +  B)  =  cos  A  cos  B  -  sin  A  sin  B 

Let  h  +  k  +  vC  =  2m      k  =  2m  =  i  -  h 

Let  us  look  at  h  +  k  +  i  =  2m.  The  magnitude  of  m  is  unimportant;  however, 

whether  it  is  odd  or  even  is  important.  Three  conditions  can  hold: 

If  h  is  even,  X  is  odd  and  k  is  odd.  ■ 
If  h  is  even,  X  is  even  and  k  is  even. 
If  h  is  odd,  J.  is  odd  and  k  is  even. 

But  we  notice  that  this  is  the  same  behavior  as  exhibited  by  h  +  >£,  so 

2m  -  h  -xmay  be  replaced  by  h  +X  everywhere  except  in  the  cosine  terms.  To 

deal  with  these  possible  variations,  the  cosine  terms  were  expanded 


so 


^M"^^  i^lS  U,X)     w-O*  ^^1^1 


From  this  point  it  is  a  matter  of  combining  the  conditions 


h  +  k  +  1  = 


i 


ceven 
even   h  +  k  +  l  =  2n,  m  (_odd 

odd   fh  +  k  +  1  =  S  =  (1,  ^,  9,  13 )  =  u 

\h  +  k  +  1  =  S  =  (3,  7,  11,  )  =  V 

even   h  =  2P,  P 


Ceven 
|odd 

odd 


^■"^"^    toid 


even       y «  y«  ^     \  even 


odd 


hO 


DETERMINATION  OF  AlCLES  BETWEEN  PLAIffiS 

Calcite  is  rhombehedral  in  structure.  This  implies  that  the  three 
periodic  distances  are  equal,  i.e.,  aj_  =  a2  =  a,,  and  the  angles  betv;een  the 
cell  axes  are  the  same.  For  calcite  the  angle  is  U6°07 '  and  a,  the  periodic 
distance,  is  6. 361  A°.  When  dealing  with  a  crystal,  it  is  easier  to  examine 
the  reciprocal  lattice  rather  than  the  natural  lattice.  In  the  reciprocal 
lattice,  each  plane  of  the  crystal  is  represented  by  a  vector  H,  which  is 
perpendicular  to  the  plane  and  has  magnitude  equal  to  -r  for  that  plane.  A 
periodic  distance  in  reciprocal  lattice  is  defined  as 

a.  X  a, 


(0         \-? 


a. -aj  X  a^ 


Since  the  natural  lattice  is  not  cubic,  tetragonal  or  orthorhombic  neither 
is  the  reciprocal  lattice.  Consequently,  the  primary  axes  are  not  orthogonal 
and  ^1*^2*  1^1 1  r?l  *  "'""  order  to  calculate  angles  between  planes  it  was 
first  necessary  to  find  the  angle  between  the  primary  axes  in  the  reciprocal 
lattice.  This  computation  was  carried  out  in  the  following  manner. 

From  geometrical  consideration  the  scalar  product  between  b2_*b2  for  a 
rhombehedral  parallelopipead  is  . 

fcOS^'<  -  COS'K) 


1   2   a^d  +  2  cos^iK  -  3cos2m) 

since  the  magnitude  of  b-,  =  bp  "  't>-5>  then  b-j_*  bg  =  b,*"  cos/jj  where  /a,   is  the 
angle  between  b-^^  and  b2. 

p         sin^-K 
^  ^    ^         a^d  +  2cos3»<  -  3cos2-\) 


/ 


U2 


so 


f\^^  cos/2    =  ,  "^"^"^    =  -0.U09 

With  this  information  the  interplanar  angles  were  determined. 

To  compute  the  angles  between  the  various  planes,  the  scalar  products 
between  the  various  planes  were  considered. 

H(hkl)  •  H(h'k'l')  =  (hbi  +  kb2  +  lb3)  •  (h'b^  +  k'b2  +  l'b3) 

=  hh'b-j^  •  b-]_  +  kk'b2  '  b2  +  ll'b3  -  b^ 
hk'bi  '  b2  +  hl'b;i_  •  b^  +  kl'b2  *  bi 
kl'b2  •  b3  +  Ih'b^  •  b^  +  Ik'b^  •  b2 

=  hh'b-L2  +  kk.'b22  +  ll'b32  +  hk' (bi2)  (_o.a09)  +  hi' {h-^)^i -0,^09) 
+  kh'(bi2)(.o.U09)   +  kl'(bi2)(_o.U09)  +  Ih' (bi2)(-O.U09) 
+  lk'(b;L^)(-0.U09) 

%kl   •  %'k'l'   =    [(^l^'   "^  ^^'   +  11')  +  fhk'   +  hi'   +  kh'   +  kl'   +  Ih'   +  Ik'] 

[-0.h09]]bi2 

but     %ki  •  Hj^.k'l'    =    |%kll    K'k'l'l        cos  ^ 

cos  ca  =     "hkl    •   Hh'k'l' 
'°'^        |Hhkl|     |Hh.k.l'|. 

|%kl|  =  f h^bjL^  +  k2b22  +  l2b32  +  2hkb3_   •  b2  +  2hlb;i^   •  b3   +  2klb2   •   h^ 

=  b^  [h2  +  k2  +  l2  +  (-0.U09)    (2hk  +  2hl  +  2kl)]     ^/^ 

cos  0  = 

[]{hh'kk'+ll']     +/"hk'+hl'+kk'+kl'+lh'+lk'}r-0.[|09]  ] 

fh2+k2+l2+(.o.U09)(2)(hk+hl+kl)]l/2^h,2+k'2+it2+(_o.ao9)(2)(h'k'+h'l'+k'l')f/2 

As  an  example,  the  angle  between  the  (211)  and  (121)  planes  is 

cos  0   =  ff2  ^  2^  l|  ^  a  ^  2  ^  1  ^  1^  1  ^  2]r(-O.U09)]]   ^  0^ 
^°^  ^  [U  +  1  +  1  -  (O.U09)(2)(2  +  2  +  1}  I  1.91 

so  0  =  75 'U^ 
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When  choosing  the  planes  to  examine,  two  criteria  were  taken  into 
consideration.  First,  it  is  desired  that  the  crystal  structure  factors  be 
different  so  that  an  attempt  could  be  made  to  isolate  contributions  of  the 
individual  atoms.  Since  the  Debye  Waller  factor  contains  a  vibration 
constant,  it  was  desired  that  the  reciprocal  lattice  vectors  chosen  be 
almost  parallel,  that  is,  the  angles  between  them  be  small.  With  these 
criteria  in  mind  and  also  considering  intensities  of  diffraction,  the  211, 
I422,  h23,  301,  101,  and  213  planes  were  chosen. 


X-RAY  DIFFRACTION  TECHNIQUES 

In  order  to  determine  experimentally  epsilons  for  the  planes  in 
calcite,  the  crystal  had  to  be  aligned  in  such  a  manner  that  the  reciprocal 
lattice  vector  for  the  desired  plane  lay  in  the  plane  of  incidence  of  the 
x-rsys;  that  is,  the  plane  defined  by  the  intersection  of  the  incident  beam 
and  the  diffracted  beam.  In  addition,  the  reciprocal  lattice  vector  must  be 
at  such  an  angle  that  it  bisects  the  angle  between  the  incident  and 
diffracted  beams.  In  order  to  do  this,  the  previously  determined  angle 
between  the  face  planes  and  each  reciprocal  lattice  vector  were  used. 

The  crystal  was  oriented  such  that  the  face  plane  121  was  perpendicular 
to  the  x-ray  beam  when  the  diffractometer  was  set  for  zero  degrees.  The 
crystal  was  rotated  fifteen  degrees  (the  compliment  of  the  angle  between 
the  121  and  the  211).  Next  the  diffractometer  was  set  on  the  Bragg  angle 
for  the  211.  Then,  without  a  slit  or  collimator  on  the  detector,  the  angle 
X  was  varied  until  the  diffraction  maximum  was  reached.  At  this  point  a 
collimator  and  slit  were  placed  on  the  detector,  and  the  angle  9.   was  varied 
to  obtain  a  maximum.  The  crystal  was  then  aligned  foi*  diffraction  from 
both  the  ko<  and  the  k/9t  wavelengths  of  the  radiation.  The  diffractometer 
was  motor  driven  to  rotate  through  all  Bragg  angles  at  a  rate  of  0.2  degrees 
per  minute.  The  areas  of  the  graphs  for  the  intensity  diffracted  by  k  *< 
and  k/3>  wavelengths  for  each  set  of  planes  were  the  data  from  which  epsilon 
was  calculated. 

The  same  slit  was  not  used  for  every  plane.  For  most,  the  0.1  degree 
slit  was  used,  but  for  the  very  weakly  diffracting  planes,  the  0.3  degree 
slit  was  used.  This  was  done  since  the  very  weak  planes  did  not  have 
sufficient  intensity  to  be  counted  with  the  0.1°  slit.  Since  the  ratio  of 
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the  areas  under  the  k»<  and  k/i.   diffracted  peaks  was  not  affected  by  the 
size  of  the  slit,  there  was  nothin^i;  loot  in  going  to  the  larger  slit. 

In  order  that  the  diffraction  graphs  be  limited  to  the  k°<    and  k/2, 
radiation  only,  the  base  line  on  the  pulse  height  analyzer  was  set  so  that 
the  k°<  line  was  just  barely  clipped.  This  insured  that  no  wavelength 
longer  than  the  k°<   was  counted.  The  window  was  closed  down  until  the  k /^ 
line  was  just  barely  clipped.  This  insured  that  no  wavelength  shorter 
than  the  k/2.   was  counted. 

In  the  data  analysis  it  was  necessary  to  measure  the  area  of  the 
intensity  curves  for  the  incident  k  "^  and  k/S.  radiations,  Iq^  and  Iq^   , 
respectively,  and  take  their  ratios.  To  obtain  this,  a  powdered  quartz 
sample  was  used.  The  Bragg  diffraction  peaks  for  the  most  intensely 
diffracting  planes  of  quartz  powder  were  used.  The  ratio  of  the  incident 
intensity  of  the  k*^  and  k^  wavelengths  (Iq^  Ao/a  ^  ^^®  related  to  the  area 
under  the  graphs  of  the  powder  quartz  sample  by  the  equation 

where  A«.  and  A^    are  the  areas  under  the  k-^  and  k/3.  Bragg  reflections, 
respectively,  0^  and  0^  are  the  Bragg  angles  for  the  respective  wavelengths 
in  quartz,  and  1^^     and  Iq^  are  the  intensity  of  the  incident  k<x   and  k/^ 
wavelengths,  respectively. 


ANALYSIS  OF  DATA 
In  the  analysis  of  data  the  followinf;  terms  were  used: 

^^  is  the  diffraction  angle  for  the  K  °<  line. 

e^  is  the  diffraction  angle  for  the  K/3  line. 

A  is  the  angle  between  the  reciprocal  lattice  vector  of  the  face 
plane  of  the  crystal  (121),  and  the  reciprocal' lattice  vector 
of  the  plane  from  which  the  diffraction  occurs. 

n.^  -  /-«,  A  ^ 
n,  --  M-  ^«  ^ 

X.«  is  the  intensity  of  the  incident  K<K  wavelength  as  defined  before. 

I„^  is  the  intensity  of  the  incident  K /a,  wavelength. 

X^(ha),I.(RaPo')  iToC*^*) ,  r«  (''^- P«<)    are  modified  Bessel  functions  for 
specific  arguments.  All  Bessel  functions  are  zero  order. 

Dreiling  and  Dragsdorf?  showed  that  the  integrated  intensity  of  an  x-ray 

diffracted  beam  is  equal  to,  for  a  single  polarization,  P  =  cos  20 

f H  -  '^  .^^(-A^.O  ro(p.^O   cj3aiiJil 


For  the  case  of  two  polarizations  being  present  in  the  incident  beam,  the 
form  for  the  integrated  intensity  is  more  complicated.  There  are  na^  two 
terms,  each  of  which  contains  a  Bessel  function  of  a  different  argument. 
The  value  of  epsilon  can  be  obtained  by  an  iteration  procedure.  In  the 


.,*>'•■■  ?^-/ 
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equation, 

A—  a.  ^0  L 
I,  (  H.  '''  «^  t^  !l  e-0  )  J 

f-,  and  fp  are  the  fractional  parts  of  each  polarization  in  the  incident  beam. 
For  this  experiment  and  unpolarized  radiation  ^1  ="  ^2  '°  '^1'^'     -^^  '^''^®  ratios 
of  the  integrated  intensities  of  the  K<^  and  KA  line  are  considered,  the 
following  equation  results: 


a: 


j^Ji^r\,   ^  ^"a 


ja^~~r.    X 


>!££. 


JhiP. 


(-I"l.) 


1 


AyjLp 


(-f^^) 


All  of  the  above  parameters  except  e  are  known.  For  the  solution  the  equation 


was  rewritten  in  the  form 


ras  rewritten  m  1 

(b)   l^f- 


-  [Io(^k)+  ''n'rofp,  K.)] 


This  equation  originally  was  solved  graphically  (see  figure  16).  A  graph 
consisting  of  curves  for  both  the  left  and  right  terms  of  the  equation  was 
drawn  and  the  intersection  gave  the  sought  value  for  epsilon.  However,  the 
method  was  time  consuming  so  a  computer  program  was  written  to  solve  for 
epsilon. 

In  order  to  solve  the  equation  by  use  of  a  computer,  an  iterative 
process  is  used  in  which  a  value  of  e  is  picked  between  0.0  and  1.0.  If 
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the  left  side  of  the  equation  is  larger  than  the  right  side,  the  value  of 
epsilon  is  increased  by  0.1  and  the  terms  are  evaluated  again.  The  process 
is  repeated  until  the  right  side  becomes  larger  than  the  left,  at  this 
time  the  value  of  epsilon  is  decreased  by  0.01.  This  process  is  continued 
until  a  final  value  of  epsilon,  accurate  to  0.001,  is  obtained.  A  listing 
of  the  program  is  in  the  appendix. 

Table  U  shows  the  variation  in  the  resulting  value  of  e  as  computer 
calculated  for  various  assumed  variations  in  the  parameters  of  the  equation 
relating  the  K°<  to  the  K/^  diffracted  intensities.  That  is,  equation  13  is 
rewritten  here. 


^ ^ 


'^=  ^    JT^^-T?  ^^H?-  -^(^•"a^^-Z^-^O-^^CAi.-/^.) 


has  variables  with  potential  errors.  The  table  thus  indicates  directly  the 
variation  of  e  for  the  fractional  variation  in  the  individual  parameters. 

The  measurements  of  t^  were  taken  five  times  with  an  average  derivation 
of  0.0005  i^™  0^  0'^  per  cent.  Two  graphs  each  of  the  diffraction  peaks  were 
taken.  In  addition  each  peak  was  measured  a  minimum  of  three  times.  The 
largest  percentage  difference  between  any  two  measurements  of  the  same  line 
was  8  per  cent  for  the  (U22)  plane.  For  the  other  planes  the  percentage 
difference  was  less  than  3  per  cent.  The  values  for  jj^  and  Vq  were  taken 
from  Cullity^  and  may  be  considered  accurate  to  U  per  cent. 
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TABLE  U 


ERROR  ANALYSIS 

Parameter 

Change  in  Parameter 
(Percentage) 

Change  in  Epsilon 
(Percentage) 

P 

together 

5^. 
6  % 

5^ 
2.57. 
2.8  ^ 

CONCLUSION 

A  comparison  between  the  values  of  epsilon  for  the  two  theoretical 
and  the  experimental  values  of  epsilon  can  now  be  made.  The  values  for 
these  three  are  presented  in  table  $. 

TABLE  $ 

THEORETICAL  AND  EXPERIMENTAL  VALUES  OF  EPSILON 


Plane 

Calculated 
Formal 

Calculated 

H-r 

Experimental 

211 

0.961 

0.9U2 

0.872 

310 

0.919 

0.883 

h22 

0.919 

0.883 

0.662 

lOl 

0.931 

0.883 

0.807 

213 

0.93h 

0.929 

0.866 

U23 

0.019^ 

0.108 

Several  items  are  worthy  of  note.  First,  the  experimental  values 
for  epsilon  are  all  lower  than  the  calculated  values.  The  stronger  dif- 
fracting planes  are  farther  from  agreement  than  the  weaker  diffracting 
planes  (the  U22  is  an  exception).  This  was  not  particularly  startling, 
for  exactly  the  same  thing  occurred  in  the  measurements  that  Mwanje^'^  made 
on  quartz.  In  table  6  the  ratio  of  the  experimental  value  of  epsilon  to 
the  formal  theoretical  value  is  given  for  each  plane. 


i^A 
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.     TABLE  6 

THE  RATIO  OF   CALCULATF.D  EXPtlRIMENTAL 

VALUES  OF  EPSILON  AND  ROOT  MEAN 

SQUARE  DISPLAGEI^JiNTS 

Plane  e^/e  <R*//^ 

211  1.10  0.37 

U22  1.3b  0.3)a 

101  1.16  0.36 

213  1.0?  0.20 

Since  the  experimental  and  the  calculated  values  of  epsilon  are 
related  by  the  equation  . . 

3        %       J  ) 

this  ratio  is  simply  a  measure  of  the  Debye  Waller  Factor  M.  Since  this 
factor  is  different  for  the  different  planes,  herein  lies  an  explanation 
for  the  low  values  of  epsilon. 

The  two  theoretical  calculations  should  be  different,  because  they 
represent  two  different  ideas  about  the  contribution  of  the  various  atoms 
to  the  total  absorption  of  x-rays  in  the  crystal.  Actually  the  H*r  method 
should  give  the  highest  values  for  epsilon,  because  it  considers  the  atom 
2S  a  point  absorber.  The  formal  theory  considers  the  absorption  by  weighing 
each  electron  shell  of  each  atom  separately,  and  in  doing  so,  produces 
lower  values  for  epsilon.  The  fact  that  for  this  crystal  this  does  not 
seem  to  hold,  is  due  to  a  difference  in  the  calculation  for  the  total 
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absorption  coefficient  for  the  crystal.  The  H'r  rr.ethod  yields  a  value  of 
Pq   =  205/cin.  for  copper  K<^  radiation,  whereas  in  the  formal  theory  the 
vslue  for  /^^  is  196/cm.  If  a  common  value  for  jL>  q  ^^  used  for  both 
calculations,  then  the  H'r  method  will  yield  the  higher  values  as  predicted. 

One  other  thing  that  could  cause  the  experimental  values  for  epsilon 
to  be  low  is  that  the  crystal  may  not  have  been  absolutely  perfect.  It 
certainly  did  not  yield  as  perfect  a  topograph  as  did  the  silicon  crystals 
that  Dreiling5  worked  vjith  nor  as  good  as  the  quartz  that  Kwanje^O  used. 

That  the  value  for  epsilon  is  very  dependent  on  crystal  perfection  is 
borne  out  by  a  set  of  similar  experiments  done  on  another  crystal  that  had 
a  noticeably  imperfect  topograph.  The  values  for  epsilon  were  almost 
70  per  cent  lower  than  the  calculated  values. 

It  would  seem,  then,  that  the  Lang  method  of  determining  crystal 
perfection  is  adequate  to  insure  the  establishment  of  the  wave  fields  in 
the  crystal. 


SUFiMARY 

Anomalously  large  transmitted  and  diffracted  x-ray  beams  hrive  been 
observed  as  propagated  through  thick  perfect  crystals.  These  beams  are 
anomalously  large  since  the  absorption  does  not  follov/  the  ordinary 
expression  I  =  l^t^l-tJ^t)     ,  but  more  closely  obeys  an  expression  of  the  type 
I  =  1  ■^p\^y'*'^^  ~   ^{Jwhere  epsilon,  e,  is  a  parameter  that  is  obtained  from 
the  Dynamical  Theory  of  X-Ray  Diffraction.  Epsilon  varies  from  zero  to  one 
depending  on  the  particular  type  of  crystal  and  from  which  plane  in  that 
crystal  the  diffraction  takes  place. 

Vibration  of  atoms  along  a  particular  direction  in  the  crystal  is  also 
related  to  the  value  of  epsilon  by  the  equation 

where  e  is  the  experimental  and  Eq  is  the  theoretical  value  of  epsilon. 
The  mean  square  amplitude  of  lattice  vibration  is  given  by  ,/i^  />  ,  while 
e„  is  the  Bragg  angle  and  X  is  the  wavelength  of  incident  radiation. 

Naturally  occurring  crystals  of  calcite  were  used  for  a  determination  of 
epsilon.  Crystal  structure  perfection  was  determined  by  use  of  the  Lang 
Topograph  Technique.  The  crystals  were  polished  in  order  to  remove  surface 
imperfections . 

The  experimental  value  of  epsilon  for  a  given  plane  was  determined  by 
using  the  ratio  of  the  diffracted  K  «^  to  the  diffracted  K/9  integrated 
intensities.  The  K -K  and  K^  integrated  intensities  were  related  by  an 
equation  involving  a  sum  of  two  modified  Bessel  functions  of  zero  order. 
A  program  was  written  to  obtain  an  iterative  solution  to  the  Bessel  equation. 
The  values  of  epsilon  for  the  211,  U22,  213,  and  101  planes  obtained  in 


this  manner  were  O.872,  0.662,  0.866,  and  0.80? ,  respectively. 

A  comparison  of  the  experimental  values  to  the  theoretical  value  of 
epsilon  yielded  the  root  mean  square  displacement  of  the  atoms  alon?^  the 
normals  to  the  101,  U22,  213,  and  211  planes.  The  values  of  (u'^/ 
obtained  were  0.36,  0.3U,  0.20,  and  0.37,  respectively. 
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COMPUTER  PROGRAM  FOR  CALCULATION  OF  EPSILON 

The  following  program  was  written  and  used  for  the  calculation  of 

epsilon  from  experimental  data.  An  iterative  process  was  used.  An  equation 

with  epsilon  dependent  terms  on  both  sides  was  solved  by  assiiming  an  initial 

value  of  epsilon  and  incrementing  this  value  until  epsilon  was  found 

within  0.001. 

On  the  data  card  the  following  values  were  read  in 

UA     The  value  of  the  linear  absorption  coefficient  for  Kc\   radiation 

UB     The  value  of  the  linear  absorption  coefficient  for  K/^  radiation 

TO     The  thickness  of  the  crystal  measured  with  a  micrometer  caliper 

WL     The  ratio  of  the  K«^  to  K/X  wavelengths 

OBA    The  ratio  of  the  incident  K/S  to  Ko<  x-ray  intensities 

A     The  angle  between  the  normal  to  the  surface  (that  the  x-rays 
are  incident  upon)  and  the  reciprocal  lattice  vector  of  the 
diffraction  plane 

0A2    Twice  the  Bragg  angles  in  degrees  for  the  K  «^  and  K  (^  wavelengths, 
0B2    respectively 

EA     The  areas  under  the  K°<  and  KA  diffraction  peaks,  respectively 

EB 

IPL    The  numerical  value  of  the  reciprocal  lattice  vector  for  the 
diffraction  plane 

CRYST  The  crystal  identification 

RAD    The  type  of  target  in  the  x-ray  tube 
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PROGRAM  TO  COMPUTE  EPSILON 


RDD,RUN=FREE,TIME=10,PAGES=30,LINES=60,KP=29 

00005  FORMAT ( 10 (F7.3)) 
00010  FORKAT(hOX,F9.6) 
00015  F0Rr4AT('    ',    '    TIffi  VALUE  OF  EPSILON  OBTAINED  FORM  ANO^:ALOUS  TRAKSMI 

ASSIGN  OF  X-RAYS  FOR  THE   ',2Ah,'    PLANES  OF    • ,UAU)    ' 
00017   FORKATC    ','  USING    ',A2,'    RADIATION  IS    ',F8.6) 

DH'iENSION  IPL(2),   CRYST(U),   RAD(1) 
00020  F0Rl'iAT(2AU,hAh,A2) 
00025  F0RMAT(10X,U(F15.6)) 
UOOl  FORMAT ( ' 1 '  ,  18 ( AU ) )    ' 
UOl  FORMAT('0',UX,3(F15.6)) 
00035  FORHATCO',2X,'UA',5X,'UB',5X,'TO',5X,'WLS5X,'OBA',5^,'A<,5X,'OA' 

*,5-X, '  OB'  ,6X,  'EA>  ,5X,  'EB' ) 
UOO  F0imAT(l8(AU)) 

AA=0.25- 

BB=A A/9.0 
C=BB/16.0 
D=C/25'.0 
E=D/36.0 
FI=E/U9.0 
S=FlM.O 
T=S/81.0 
U=T/100.0 
V=U/121.0 
W=V/lUh.0 
YY=W.Vl69.0 

db'Iension  title (18) 

210  READ (1, UOO) title 

READ(1,5)UA,UB,T0,WL,0BA,A,0A2,0B2,EA,EB 

pj;ad(i,20)ipl,cryst,rad 

0A=0A2/(2-.0^57-295) 
0B=0B2/(2.0*57.295) 

a=a/57.295 

TA=(T0/2.0)*(((l-.0/SIN(0A+A))-(r.0/SIN(0A-A)))) 
TB=(TO/2.0)^'^(((1.0/SIN(OB+A))-(1.0/SIN(OB-A)))) 

UATA=UA*TA 
UBTB=UB-;^TB 
CA=C0S(0A) 
CB=C0S(0B) 
IF(UA)212,213,212 
212  PA=2*CA-i'<-;^2-1.0 
PB=2^:-CB*-^'-2-1.0 
AP=UATA^PA 
BP=UBTB^PB 
ARU=(UATA-UBTB) 
EX=EXP(ARU) 
CONST=(EA4fOBA^iWL-«-CA*EX)/(EB*CB) 

bbb=sin(ob+a)/sik(a-ob) 
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PROGRAM  TO  COMPUTE  EPSILON   (continued) 


BAA=SIN(OA+A)/SIN(A-OA) 

C0NSTs=C0NST^(-(BBB/BAA) 

HIO=0.0 

EP=0.5 

1=0 

L=0-  ■       ■ 

F=0.1 

FF=0.01 

FFF=0.001 

OA=OA-^2.0-"-57.295 

0B=0B-s-2. 0*57. 295 

WRITE (^,h00l) TITLE 

WRITE (3,^5) 
WRITE(^,5)UA,UB,T0,WL,0BA,A,0A,0B,EA,EB 

DO  50  J=1,?0 

WRITE (^, 10 )EP 

W  =  (UATA-"-EP)^;-«-2/U.0 

X=(AP-"-EP)^:~-2A.O- 

Y=(UBTBitEP)-"ri;-2/U.O 

Z=(BP-5<-EP)-)Hc-2/h.O 

CI0=1.0+W-K-(l.0+W*(AA+W^:-(BB+W^t(C+W*(D+W*(E+W*(FI+W*(S+W*(T+V7*(U+W*( 
W+W^KVA^+W-x-YY)  ))•))))))))) 

DI0=PA(1.0+X-;:-(l.0+X*(AA+X-5^(BB+X^t(C+X*(D+X*(E+X*(FH-X*(S+X*(T+X-^-(U 
-;^+X (V+X^KW+X^cYY) )))))))))))) 

EI0=1.0+Y-x-(l.0+y*(AA+Y^i-(BB+YiKC+Y^'i-(D+Y-^-(E+Yrr(FI+Y^<-(S+Y*(T+H^(U+Y*( 
-x-V+Y^KVM+Y^-YY)  )))•)))))))) 

FIO=PB^Kl.O+Z^r(l.O+Zrr(AA+Z^(BB+Z*(C+Z*(D+Z*(E+Z(FI+Z*(S+Z*(T+Z*(U 
r<-+Z^KV+Z-!>KW¥+Z-i^YY) )))))))))))) 

DI0=DI0-}f-PA-iH^2 

FI0=FI0-;;-PB^h;-2 

WRITE(?,25)CIO,DIO,EIO,FIO 

CIO=CIO+DIO 

EIO=CONST^;-(EIO+FIO)  '     ' 

GI0=CI0-EIO 

WRITE  (■=!  ,U01)  CIO, EIO,GIO 

IF(L)190,58,190 
00058  IF(I)70, 60,70 

60  IF(HI0)150,100,130 
100  IF(GI0)120,200,110 
110  EP=EP-F 

HI0=GI0 

GO  TO  50 
120  EP=EP+F 

HIO=GIO 

GO  TO  50 
130  IF(GI0)1U0,200,110 
lUO  EP=EP+FF 

1=1 
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PROGRW  TO  COMPUTE  EPSILON   (continued) 


HIO=GIO 

GO  TO  50 
150  IF(GI0)120,200,160 
160  EP=EP-FF 

1=1 

HIO=GIO 

GO  TO  50 
70  HR=GIO/HIO 

IF (HIO) 17 0,200, 180 
170  IF(HR.)175,200,1U0 
175  EP=EP-FFF 
.    L=l 

HI0=GI0 

GO  TO  50 
180  IF(HR)185,200,160 
185  EP=EP+FFF 

L=l 

HI0=GI0 

GO  TO   50 
190  HR=GI0/HI0 

IF(HI0)192,200,19L 
192   IF  (HR-) 200, 200, 185 

19U  IF (HR) 200, 200, 175 
200  \\^ITE(?,10)EP 

¥PITE(3,15)IPL,CRYST 

VJRITE(  1,17)  RAD,  EP 

GO  TO  210 
50  CONTINUE 

WRITE(3,10)EP 

WRITE  ( -^,10  )EP 

WRITE(3,15)IPL,CRYST 

WRITE (3, 17)  RAD, EP 

GO  TO  210 
213  STOP 

END 
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Anomalously  large  transmitted  and  diffracted  x-ray  beams  have  been 
observed  as  propagated  through  thick  perfect  crystals.  These  beams  are 
anomalously  large  since  the  absorption  does  not  follow  the  ordinary  expression 
I  =  Iq  -f'o't,  but  more  closely  obeys  an  expression  of  the  type 
1  =  1  -A'o't'^l  -  ^)  where  epsilon,  e,  is  a  parameter  that  is  obtained  from  the 
Dynamical  Theory  of  X-Ray  Diffraction.  Epsilon  varies  from  zero  to. one 
depending  on  the  particular  type  of  crystal  and  from  which  plane  in  that 
crystal  the  diffraction  takes  place. 

Vibration  of  atoms  along  a  particular  direction  in  the  crystal  are  also 
related  to  the  value  of  epsilon  by  the  equation 

where  e  is  the  experimental  and  Cq  is  the  theoretical  value  of  epsilon. 
The  mean  square  amplitude  of  lattice  vibration  is  given  by  ^«-l V  ,  while 
0^   is  the  Bragg  angle  and  \  is  the  wavelength  of  incident  radiation. 

Naturally  occurring  crystals  of  calcite  were  used  for  a  determination  of 
epsilon.  Crystal  structure  perfection  was  determined  by  use  of  the  Lang     .. 
Topograph  Technique.  The  crystals  were  polished  in  order  to  remove  surface 
imperfections. 

The  experimental  value  of  epsilon  for  a  given  plane  was  determined  by 
using  the  ratio  of  the  diffracted  K**'  to  the  diffracted  K/Sl  integrated 
intensities.   The  K<<    and  K/SL  integrated  intensities  were  related  by  an 
equation  involving  a  sum  of  two  modified  Bessel  functions  of  zero  order. 
.  A  program  was  written  to  obtain  an  iterative  solution  to  the  Bessel  equation. 
The  values  of  epsilon  for  the  211,  U22,  213,  and  101  planes  were  obtained  in 
this  manner.  The  values  of  epsilon  obtained  were  o.vrs'^  o,6»7^  o.iyo  ""^ 

[.-;,:■■■.  '        ■  "    '■  ■  ■         _      ■  ■ 


A  comparison  of  the  experimental  values  to  the  theoretical  value  of 
epsilon  yielded  the  root  mean  square  displacement  of  the  atoms  alon^^  the 
normals  to  the  lOT,  U22,  213,  and  211  planes.  The  values  of  u  obtained  were 


